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Abstract We present a comprehensive review of past research into adiabatic quan- 
tum computation and then propose a scalable architecture for an adi- 
abatic quantum computer that can treat NP-hard problems without 
requiring local coherent operations. Instead, computation can be per- 
formed entirely by adiabatically varying a magnètic field applied to all 
the qubits simultaneously. Local (incoherent) operations are needed 
only for: (1) switching on or off certain pairwise, nearest-neighbor in- 
ductive couplings in order to set the problem to be solved and (2) mea- 
suring some subset of the qubits in order to obtain the answer to the 
problem. 
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1. Introduction 

Adiabatic quantum computation 1 is a recently proposed, general ap- 
proach to solving NP-hard combinatorial minimization problems. 2 It 
consists of constructing a set of qubits with a time-dependent Hamil- 
tonian H (i) whose starting point H s has a ground state that is quickly 
reachable simply by cooling and whose final point H p has couplings that 
encode the cost scheme of a desired minimization problem. The name 
"adiabatic" comes from the fact that if the qubits are initialized in the 
ground state of H s and if H(t) is varied slowly enough, then the qubits 
will overwhelmingly be in the ground state of H(t) at all times t, thus 
in principle completely bypassing the usual concern about local minima 
in H p confounding the search for the problem's solution. 
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Sections 2 and 3 review the literature to date on adiabatic quantum 
computation. Section 4 presents the main result of this paper: a scalable 
architecture for an adiabatic quantum computer. Section 5 concludes. 

2. Literature Review: The Time Complexity of 
Adiabatic Quantum Computation 

The general time-complexity of adiabatic quantum computation is 
still an open problem. Intuitively, the key to fast, successful adiabatic 
computation is to ensure that H(t) always possesses a sizable spectral 
gap between its instantaneous ground and excited states, thus allowing 
one to vary H(t) quickly while still not ever providing enough energy 
to excite the system significantly out of its instantaneous ground state. 
Determining the asymptotic time complexity of adiabatic quantum com- 
putation remains an open problem because for n qubits, H(t) is a 2 n x 2 n 
matrix, and hence analytically obtaining lower bounds on the spectral 
gap is extremely dimcult unless H(t) possesses significant symmetry. 

Presently, there are three major results on the time complexity of adi- 
abatic quantum computation. Together, they yield the general picture 
that for essentially unstructured NP-hard minimization problems like 
finding the ground state of a random-field magnet, adiabatic quantum 
computing should offer polynomial speedups — probably in the range of 
cube-root to sixth-root — over a Metròpolis algorithm (simple cooling) 
approach to the same problems. Furthermore, for structured NP-hard 
minimization problems such as 3-SAT and MAX CLIQUE, adiabatic 
quantum computation can quite possibly offer exponential speedups over 
any classical algorithm (assuming P ^ NP) — at least on average. 

More specifically, the first major time-complexity result is that adia- 
batic quantum computation can search an unsorted list of N items for 
a single identifiable item in 0(s/N) time, 3 thus matching the speedup 
of Grover's algorithm. 4 This, however, is a naïve way of characteriz- 
ing the performance of adiabatic quantum computation on unstructured 
searches because an nth-root speedup still means that the minimum 
spectral gap in H (í) is exponentially shrinking, and practically one al- 
ways needs to operate above some finite temperature limit. 

It is therefore much more realistic when analyzing unstructured prob- 
lems to abandon the idea that the system can always remain overwhelm- 
ingly in the instantaneous ground state of H{t) and instead view adia- 
batic quantum computing as an enhanced form of annealing. The poten- 
tial for enhancement intuitively derives from two facts. One, quantum 
systems can tunnel through energy barriers instead of waiting for ther- 
mal excitation over them. Two, with a time dependent Hamiltonian 
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H(t), one need not face the same barriers as those in energy landscape 
H p of one's problem. 

For example, consider an Ising model of the form 

where the couplings {Jij} are randomly drawn with a zero-mean Gaus- 
sian distribution from the interval [—1,1]. If the Ising model cannot be 
drawn as a planar graph (with spins as vertices and nonzero couplings 
as edges), then finding the ground state of H p is an NP-hard problem. 5 
Quite general considerations of the distribution of low-lying excitations 
of any frustrated system such as that of Equ. (1) imply that simply cool- 
ing the system from T = ootoT = 0ina time r (physically, "oo" means 
"sufficiently high to make all possible states essentially equiprobable" ) 
leaves an average residual energy 

(E(t)) - E ground ~ ln(r)-* (£ < 2) (2) 

in the limit of large r. 6 Monte Cario simulations imply that £ can be as 
low as 1 for cooling to the ground states of some NP-hard random-field 
Ising models. 7 In contrast, if instead one applies a transverse magnètic 
field to the system that decreases from H = oo to H = in time r (all 
done at T = 0), then according to a model 8 that calculates the residual 
energy as the cumulative result of Landau-Zener type transitions 9 at the 
avoided level crossings between the instantaneous ground state and first 
excited states of the resulting time-dependent Hamiltonian 

H(t) = H P + r(t) J2 4 i} (IXO) = oo, r(r) = 0) (3) 

i 

one finds much improved performance 

(E(t)) - E ground ~ ln(r)-^ (£ Q[/ « 6) (4) 

The greater size of the exponent Çqjj in Equ. (4) over the exponent £ in 
Equ. (2) is the basis the second major time-complexity result — adiabatic 
quantum computation on a given unstructured NP problem should at 
least provide cube-root and perhaps as much as sixth-root speedups over 
a Metròpolis algorithm (simple cooling) approach to the same problem. 

The third, and perhaps most interesting, time-complexity result — 
the possibility of exponential speedup on almost all difficult instances 
of NP-hard problems — comes from explicit numerical integrations of 
Schròdinger's equation to simulate adiabatic computation on sets of 
small instances of these problems (< 26 bits is roughly the limit of 
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current supercomputers). Since the class NP is a worst-case complexity 
measure, it is important to consider problem instances that are most 
likely to be truly difficult. For example, Hogg 10 studied random 3-SAT, 
which is the problem of finding a string of n bits that completely satisfies 
m Boolean clauses, each of the form a OR b OR c where 

1 the literals a, b, and c refer to 3 distinct bits chosen uniformly at 
random 

2 each literal, with 50% probability, contains a negation. 

In the limit of large n, random 3-SAT instances exhibit a first-order 
phase transition from being highly likely to have a satisfying assign- 
ment when m/n < 4.25 to being highly unlikely to have a satisfying 
assignment when m/n > 4.25, and it has been observed that the most 
difficult instances to solve appear to fall on the phase boundary. 11 Hogg 
studied instances generated as near as possible to this phase bound- 
ary m/n = 4.25 and found that for n < 26 bits (the largest number 
that could be numerically analyzed) the median run time scaled only as 
0(n 3 ). 

Random 3-SAT instances at the phase boundary still tend to have 
exponentially many satisfying solutions which, although they constitute 
an exponentially small fraction of the 2 n possible bit strings, one might 
worry constitute enough degeneracy to distort study of such small in- 
stances of < 26 bits. Hence, the earlier studies of Farhi et al. 1 and Childs 
et al. 12 focused on randomly gener ating NP-complete problem instances 
with unique satisfying assignments ( "USA instances" ) rather than some 
ratio of constraints to bits. Median run times scaling only as 0(n 2 ) 
were observed for both USA instances of random EXACT COVER up 
to 20 bits 1 and random MAX CLIQUE up to 18 bits. 12 (Random EX- 
ACT COVER is the problem of finding a string of n bits that satisfies m 
clauses of the form a + b + c = 1 where the literals a, b, c are distinct bits 
chosen uniformly at random. Random MAX CLIQUE is the problem of 
finding the largest subgraph for which every pair of points is connected 
by an edge within an n-vertex graph constructed such that each pair of 
vertices is connected by an edge with 50% probability.) 

The above MAX CLIQUE result may be seen to carry extra signifi- 
cance because it is widely conjectured there does not exist any classical 
polynomial-time algorithm that can identify cliques in such n-vertex 
random graphs of size > (1 + e) log 2 (n) for any e > 0. In contrast, the 
maximum clique of an n-vertex random graph is almost always > log 2 (n) 
and on average is 2 log 2 (n) . This discrepancy is in fact used as the basis 
for a cryptographic scheme. 13 
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3. 



Literature Review: Robustness of Adiabatic 



Quantum Computation versus Decoherence 

Beyond the potential for significant speedups over classical comput- 
ing, which is an attraction that adiabatic quantum computing shares 
with other proposed methods of quantum computing, adiabatic quantum 
computing is especially attractive for being intrinsically robust against 
environmental noise. 14 Firstly, since an adiabatic algorithm aims to keep 
its qubits in the instantaneous ground state of the computational Hamil- 
tonian H(t), dissipation is not intrinsically harmful, but potentially help- 
ful. (Of course, one cannot depend wholly on dissipation because then 
classical thermal annealing could have solved the problem with equal ef- 
ficiency.) Secondly, if indeed the spectral gap of H(t) is large enough so 
that the system can be kept overwhelmingly in just one energy eigenstate 
(namely, the ground state), then the dephasing problem can be essen- 
tially circumvented. This is because so long as H (t) dominates over the 
environmental couplings to the qubits, the energy eigenstates of H(t) 
will be the preferred basis for dephasing, meaning that an adiabatic 
quantum computer based on H(t) would never present the environment 
with a superposition of states that it could dephase. 

Such dominance is easy to achieve in the canonical spin-boson model 
with an Ohmic environment having an exponential cutoff after a fre- 
quency w c . 15 That is, if A G denotes the bare energy splitting in the 
qubits due to coherent tunnelling between their limiting, classical states 
in the absence of coupling to the environment and a represents the (di- 
mensionless) slope of the environmental spectral density at oj = 0, then 
one can conclude non-perturbatively that the environmentally corrected 
tunnel splitting A is 



Hence, if a <C 1 and uj c is only a few orders of magnitude above A , 
then A « A D meaning that the qubits' Hamiltonian does indeed dom- 
inate over the environmental couplings and its energy eigenbasis is the 
preferred basis for dephasing. 

This desirable situation should not be changed by the addition of 
1// or other low-frequency-dominated noise to the Ohmic environment. 
Components in the environmental spectral density that are of much lower 
frequency than the bare tunnel splitting A D denote environmental dy- 
namics that are much slower than that of the qubits. Therefore, essen- 
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tially regardless of the qubit-environment coupling, these components 
can only become entangled to time-averaged properties of the qubits 
such as their energy. 16 While low-frequency components could thus sig- 
nificantly change dephasing rates, they should not change the fact that 
the environment respects the energy eigenbasis of the qubits' Hamilto- 
nian as the preferred basis for dephasing, meaning that adiabatic algo- 
rithms would still essentially circumvent the dephasing problem. 

4. Main Result: A Scalable Architecture for 
Adiabatic Quantum Computing 

Another special advantage of adiabatic quantum computing is the 
main point of this paper: there exists a simple, scalable architecture 
for adiabatic quantum computation that can handle NP-hard problems 
without requiring local coherent manipulations of the qubits. The start- 
ing point for the architecture is the fact that it is NP-hard to calcu- 
late the ground state of a 2D antiferromagnetically coupled Ising model 
placed in a magnètic field even if each Ising spin is coupled to no more 
than 3 others. 5 In other words, if G = (V, E) is a planar graph with 
vertices V and edges E and even if every vèrtex in V belongs to < 3 
edges in E, then it is NP-hard to find the ground state of 



which is an Ising model with spins at the vertices V, equal strength anti- 
ferromagnetic couplings along the edges E, and a homogeneously applied 
field giving each spin a Zeeman splitting equal to the antiferromagnetic 
coupling strength. (The direction of this applied field sets the "z-axis.") 
This proof of NP-hardness isomorphically maps the above Ising model 
onto the MAX INDEPENDENT SET problem, which is the problem of 
finding for a graph G = (V, E) the largest subset S of the vertices V 
such that no two members of S are joined by an edge from E. Solv- 
ing MAX INDEPENDENT SET for a given graph G = (V, E) is com- 
pletely equivalent to solving MAX CLIQUE for the graph's complement 
G c = (V, (V x V)/E). (In other words, G c connects with edges all the 
pairs of vertices that were unconnected in G and disconnects all the pairs 
were connected.) Therefore, the MAX CLIQUE results 12,13 referenced 
above imply that MAX INDEPENDENT SET is an equally interest- 
ing problem for treatment by adiabatic quantum computation, likely 
admitting not even a rough approximation in polynomial time classi- 
cally, yet apparently admitting efficient solution by quantum adiabatic 
evolution — at least in small instances (and thus, hopefully, for almost 
all, if not all, larger instances) . 
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There exists a simple, regular architecture that can be programmed 
to pose the above NP-hard Ising model ground state problem for any 
instance with < n spins. One way to see this is to recali that any planar 
graph G with < 3 edges per vèrtex admits a topologically equivalent 
embedding Q in a square grid of < 9n 2 vertices. 17 In other words, there 
exists a map Q that (1) maps each vèrtex v of G onto a distinct point 
Q(v) of the grid and (2) inserts "dummy vertices" so as to map an edge 
e connecting vertices u and v of G into a path Q(e) on the square grid 
such that Q{u) and Q(v) are its endpoints and such that if e and f are 
two distinct edges of G, then Q(e) and Q(í) do not intersect (except, 
obviously, if e and f share an endpoint, in which case Q(e) and Q(f) 
will also share an endpoint). 

Thus, one simple, regular architecture would be of the type depicted 
below in Fig. 1: two stacked square grids of ferromagnetic and antifer- 
romagnetic couplings, respectively, having at the center of each square 
a qubit with a ferromagnetic coupling to, say, the inductive couplings 




Figure 1. (A) A regular 4-vertex graph for an Ising model. Vertices de- 
note location of qubits, and edges denote antiferromagnetic couplings. (B) 
Topologically equivalent square grid embedding. (C) Possible implemen- 
tation. Boldface lines denote couplings switched on, dashed lines denote 
couplings switched off. 
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meeting at lower right còrner of the square. Since a useful adiabatic 
scheme would be one like Equ. (3) which merely adds a homogeneous, 
steadily decreasing transverse magnètic field to the problem Hamilto- 
nian, the coherent evolution of the qubits could be affected simply with 
a globally applied magnètic field. This avoids the need both for any 
local coherent manipulations as well as for tunable rather than merely 
switchable couplings. Furthermore, as the the switching of the couplings 
is necessary only for programming the desired instance to be solved, it 
can be done "offhne," before any coherent manipulation of the qubits 
begins, thus obviating the need for "quiet" switches. 

A final consideration for scalability is to calculate how much impre- 
cision can exist in the parameters of the qubits, the couplings, and the 
applied field before one loses the ability to encode reliably in the above 
architecture any desired n-bit instance of MAX INDEPENDENT SET. 
To answer this, first recali that Ref. 5 proved that the maximum inde- 
pendent sets of a graph G = (V, E) are in 1-to-l correspondence with the 
ground state(s) of the Hamiltonian in Equ. (6), i.e., an Ising model with 
spins at the vertices V, antiferromagnetic couplings all of equal strength 
along the edges E, and a homogeneous applied magnètic field giving 
every spin a Zeeman splitting equal to the antiferromagnetic coupling 
strength. However, so long as one is satisfied with obtaining any one of 
the maximum independent sets of G (which is reasonable since obtain- 
ing any one is NP-hard), one does not need to calibrate so precisely the 
couplings among themselves and to the applied field because the ground 
state(s) of any Ising model from the family 



where {Jj} are arbítrary positive scalars will also be ground state(s) 
of the Equ. (6) Ising model and hence maximum independent set(s) of 
G. Hence, in regard to programming the computer so that its final 
ground state has a large overlap with at least one solution of a desired 
instance of MAX INDEPENDENT SET, the proposed architecture is ro- 
bust against arbítrary variance in the magnètic moments of the qubits. 
(Note, however, that while variance in {Jj} poses no problem in pro- 
gramming the computer, it might still cause a problem in ensuring that 
the computer produces accurate solutions quickly. Specifically, as the 
variance becomes larger, the problem Hamiltonian (7) looks more and 
more like random field magnet than a "structured" NP problem, likely 
lessening the chances that the adiabatic scheme H(t) will keep a large 
spectral gap and thus permit a running time exponentially faster than 
any classical algorithm.) 
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More general errors can come from inhomogeneities in both the mag- 
nitude and the direction of the applied magnètic field and from inductive 
couplings beyond the desired ones between nearest neighbors. Thus, to 
Equ. (7) one can add the general perturbation 

h'= E EW + E E*2'W ) (8) 

k£{x,y,z} i<=V k,l£{x,y,z} %>j 



Note that Pauli operator terms such as à k or & k à\ merely prescribe 
bit and/or phase flips in the eigenbasis of the Equ. (7) problem Hamil- 
tonian, meaning that each of the perturbing terms in H' couples an 
unperturbed ground state to only one other unperturbed eigenstate. 
Therefore, if n denotes the number of qubits, A denotes the largest of 
the scalars {ej?, ò k p} in Equ. (8), and A denotes the minimum increase 
in energy caused by flipping 1 or 2 qubits of the ground state (thus, 
A ps minj[2Jj]), then Xn 2 /A serves as a simple upper bound on the 
magnitude of the flrst order correction to the ground state wavefunction. 
Higher order corrections should be well-behaved since performing 1 or 2 
bit flips on any excited state also typically entails an energy change of at 
least A, implying that the entire perturbation series correction should 
essentially converge like the geomètric series ^2 k (Xn 2 / A) k . Hence, the 
architecture's tolerance for general imprecision in the qubit couplings 
and inhomogeneity in the applied field decreases quadratically in the 
number of qubits. (If unwanted inductive couplings can be limited to 
a constant number of "next-nearest" neighbors, then the tolerance shall 
decrease only linear ly.) 

5. Conclusion 

We have reviewed the theory of adiabatic quantum computation re- 
garding its potential speedup over classical computation and its greater 
robustness toward noise over other forms of quantum computation. We 
then proposed a scalable architecture for adiabatic quantum compu- 
tation based on the NP-hardness of calculating the ground state of 
a planar antiferromagnetically coupled Ising model placed in a mag- 
nètic field. Clearly, such an architecture could be used with any qubit 
whose states have opposite magnètic dipole moments. Future work will 
present a detailed proposal for its implementation using superconduct- 
ing persistent-current qubits, 18 which constitute a current, promising 
approach to lithographable solid-state qubits. 



10 



Acknowledgments 

This work is supported in part by ARDA and DoD under the AFOSR 
DURINT Program. WMK gratefully acknowledges fellowship support 
from the Fannie and John Hertz Foundation. 

References 

1. E. Farhi et al., A Quantum Adiabatic Algorithm Applied to Random Instances of 
an NP-Complete Problem, Science 292, 472-476 (2001). 

2. e.g., M.G. Garey and D.S. Johnson, Computers and Intractability (W.H. Freeman 
and Company, New York, 1979). 

3. J. Roland and N.J. Cerf, Quantum Search by Local Adiabatic Evolution, http:/ /xxx. 
lanl.gov/abs/quant-ph/0107015 (2001). 

4. L.K. Grover, Quantum Mechanics Helps in Searching for a Needle in a Haystack, 
Phys. Rev. Lett. 79(2), 325-328 (1997). 

5. F. Barahona, On the Computational Complexity of Ising Spin Glass Models, J. 
Phys. A 15, 3241 (1982). 

6. D.A. Huse and D.S. Fisher, Residual Energies after Slow Cooling of Disordered 
Systems, Phys. Rev. Lett. 57(17), 2203-2206 (1986). 

7. G.S. Grest, CM. Soukoulis, and K. Levin, Cooling-Rate Dependence for the Spin- 
Glass Ground-State Energy: Implications for Optimization by Simulated Anneal- 
ing, Phys. Rev. Lett. 56(11), 1148 (1986). 

8. G.E. Santoro, R. Martohàk, E. Tosatti, and R. Car, Theory of Quantum Annealing 
of an Ising Spin Glass, Science 295, 2427-2430 (2002). 

9. L.D. Landau and E.M. Lifschitz, Quantum Mechanics, 3 rd rev. ed. (Butterworth- 
Heineman, Oxford, 1977), Sec. 53. 

10. T. Hogg, Adiabatic Quantum Computing for Random Satisfiability Problems, 
http://xxx.lanl.gov/abs/quant-ph/0206059 (2002). 

11. e.g., R. Monasson et al., Determining Computational Complexity from Charac- 
teristic 'Phase Transitions', Nature (London) 400, 133-137 (1999). 

12. A.M. Childs, E. Farhi, J. Goldstone, and S. Gutmann, Finding Cliques by Quan- 
tum Adiabatic Evolution, Quant. Info. Comp. 2(3), 181-191 (2002). 

13. A. Juels and M. Peinado, Hiding Cliques for Cryptographic Security, Proc. 9th 
Annual ACM-SIAM SODA, 678-684 (1998). 

14. A.M. Childs, E. Farhi, and J. Preskill, Robustness of Adiabatic Quantum Com- 
putation, Phys. Rev. A 65 012322 (2002). 

15. A.J. Leggett et al, Dynamics of the Dissipative Two-State System, Rev. Mod. 
Phys. 59, 1-85 (1987). 

16. J.P. Paz and W.H. Zurek, Quantum Limit of Decoherence: Environment Induced 
Superselection of Energy Eigenstates, Phys. Rev. Lett. 82(26), 5181-5185 (1999). 

17. L.G. Valiant, Universality Considerations in VLSI Circuits, LEEE Trans. Còmput. 
C30(2), 135-140 (1981). 

18. J.E. Mooij et al, Josephson Persistent-Current Qubit, Science 285, 1036-1039 
(1999). 



